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Recently, an ionized cloud of gas was discovered plunging toward the 
supermassive black hole, SgrA*, at the centre of the Milky Way. The 
cloud is being tidally disrupted along its path to closest approach at ~3100 
Schwarzschild radii from the black hole. Here, we show that the observed 
properties of this cloud of gas can naturally be produced by a proto-planetary 
disk surrounding a low-mass star, which was scattered from the observed ring 
of young stars orbiting SgrA*. As the young star approaches the black hole, 
its disk experiences both photo-evaporation and tidal disruption, producing 
a cloud. Our model implies that planets form in the Galactic centre, and 
that tidal debris from proto-planetary disks can flag low mass stars which are 
otherwise too faint to be detected. 



Introduction 

Observations of the galactic centre recently yielded a cloud of ionized gas and dust 
falling inward on a nearly radial orbit The cloud will reach pericenter in the summer of 
2013, approaching our galaxy's central black hole, SgrA*, at a distance of only 270 AU. As 
this plunge progresses, tidal gravity from the black hole will disrupt the infalling cloud, 
providing a unique probe of gas flow near SgrA*. 

*rmurray-clay@cfa. harvard.edu 
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The apocenter of the cloud's orbit, at rapo = 8400 AU = 0.04 pc from the black hole, 
coincides with the inner edge of the ring of young stars orbiting SgrA*, and the plane of the 
cloud's orbit coincides with that of the ring. ^ ^ The ring's age, estimated from its population 
of 0/WR stars, is ~4-8 Myr.3 At ages <3 Myr, most low mass stars host proto-planetary 
gas disks^ with radii of order 100 AU,'^ and ~l/5 of stars with mass 0.1 — IMq retain 
their disks at ages of ~5 Myr.EEl Given the black hole mass of Mbh = 4.3 x IO^Mq^ 
the tidal radius around a star of mass at a distance r = 0.04 pc from SgrA* is 

/3Mbh)^/^ ~ 4OAU(mVM0)i/3. A solar mass star could therefore host a stable 
disk with a radius of rfout ~ (it/3 ~ 12 AU on a roughly circular orbit near the inner edge 
of the young stellar ring. Similarly, an M-dwarf with mass m^, = 0.3Mq could host a stable 
disk having radius dout ~ 8 AU. 

We suggest that the newly discovered gas cloud'^ surrounds such a star, which was 
scattered away from its original ring orbit and is currently plunging toward the supermassive 
black hole (Figure 1). The star itself is too low mass to be observable, but the debris 
produced through the disruption of its proto-planetary disk allowed it to be detected. We 
first calculate the properties of the system at the cloud's observed location and argue that 
they match the observations — an ionized cloud ~100 AU in radius with density n ~ 3 x 10^ 
cm~^, an electron temperature of lO'^K, and a dust temperature of ~550K, trailed by 
a stream of gas. We then provide predictions regarding the evolution of the cloud as it 
approaches pericenter. Finally, we demonstrate that the probability of producing such an 
object is plausibly high, and we calculate the imphed rate of mass deposition by this process 
within the young stellar ring. 

Results 

Mass loss. Although the tidal radius for a solar-mass star at r = 0.04pc is 40AU, the 
tidal radius at the cloud's pericenter distance of = 270AU = lO^'^pc is only dt = lAU 
(m^/M©)^'''^. At the most recently observed epoch,^the cloud was approximately Gvp from 
the black hole, with a tidal radius of 6AU (m^/M©)^/^. Hence, the circumstellar disk is 
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already experiencing substantial tidal disruption. At the same time, the Galactic centre 
hosts an extreme flux of ionizing and Far Ultraviolet (FUV) photons. Proto-planetary disks 
in the ionizing environment near O stars in the Trapezium cluster are known to experience 
photoevaporation.'I^ The stars experience mass loss due to heating both by FUV and by 
Lyman limit photons.'^ The former heat the disk to ~10^K, generating outflows at the 
sound speed of ~3 km s~^, corresponding to the escape velocity at a distance ~100 AU 
{iti^^/Mq) from the star. Well within this distance, the FUV-driven outflow is diminished, 
though not entirely quenched.'^ At the ~10 AU and smaller distances of interest here and 
given the extreme ionizing environment, Lyman continuum (ionizing) photons dominate 
the outflow, generating a ~10^K ionized outflow moving at the sound speed of ~10 km s~^. 
This speed matches the escape velocity at a distance of d^sc ~ 10 AU (m^/Mo) from the 
star. Loss from smaller distances occurs at a reduced rate, but still generates a ~10 km s^^ 
outflow by the time the gas reaches desc- 

Which process — tidal stripping or photoevaporation — dominates mass loss from the 
disk? Currently, tidal stripping dominates the unbinding of mass from the star, and at large 
distances from the star tidal stripping determines the ultimate fate of the gas. However, 
the outflow properties of the observed cloud are nevertheless currently determined by 
photoevaporation. This can be understood as follows. Gas at a distance d > dt from 
its host star is accelerated by the tidal potential to a relative speed of Av as it moves 
of order its own radius away from its host star: Av ~ {GM-QY{/r'^){d/r){d/ Av), so that 
Av ~ t;circ(rf/r), with v^i,^ = {GMBu/rY^'^. For d = dt, Av ~ Wcirc(m*/MBH)^/^. At the 
cloud's current separation from SgrA*, r ^ 6 x 10^"^ pc ~ 1300 AU, fcirc ~ 1700 km s~^ 
and Av ~ 10 km s~^(m^,/MQ)-'^/^ at the tidal radius, comparable to the wind outflow rate. 
Hence, the motion of tidally decoupled gas is dominated by the tidal fleld. At earlier times 
in the star's plunge, Av was smaller, meaning that wind gas flowed out of the tidal radius 
faster than tidally disrupted gas. Currently, near the disk edge, the dynamics of previously 
ejected gas are set by properties of the wind, while on the ~100 AU scale of the cloud, tidal 
evolution overwhelms wind motions. 
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Although the current mass disruption rate of the proto-planetary disk, Mdis, is larger 
than the wind outflow rate, M^, wind gas emitted at earlier times dominates the currently 
observed cloud. In the short time that the infalling star has spent in an enhanced tidal field 
with dt < dout, the disk has only had the opportunity to expand by a fraction of its ~10 
AU size. At d = 8AU, the time since dt = d along the infalling star's orbit is At = Syr for 
m^, = 0.3Mq. Decoupled material has traveled only ~At'At ~ 4 AU further from the host 
star in that time. Figure 2 illustrates this point. We ask how far a test particle, released at 
a given disk radius, d, from the star when dt = d and moving only under the gravity of the 
black hole, will be from its current orbit. Exterior to the decoupled material, gas originally 
launched in a wind dominates. 

The tidal decoupling rate in the disk is Mdis ~ 2'nTjdtdt ~ 27rS(ifr(m^/3MBH)"'^^'^, where 
is the original surface density of the disk and f 2000 km s~^ at the current position 
in the cloud's orbit. For illustration, we choose a profile similar to the minimum-mass 
solar nebula: S = Eo((i/(io)^^ with Sq = 2 x 10^ g cm^ and do = lAU.'^ This choice 
yields a current Mdis ~ 3 x 10^'^ M0 yr~-'^(m*/M0)^/'^, where we set d equal to the 
current tidal radius. Photoevaporation, on the other hand, gives mass loss rate^^^*^ of 
M^ ~ 3 X 10~i°(rfout/10 K\Jf/'^{^i,ig/Dl^Y'^MQ yr-^ for disks with sizes dont > dcsc, where 
$j^49 is the ionizing luminosity, <l>i, in units of 10^^ s~^ of a source at distance D from the 
disk, with D^c = D / (Ipc). 

The photoevaporation mass loss can be derived as follows. At = desc, the 
escape velocity is comparable to the wind's sound speed Cg = 10km s~^. For d > desc, 

^ And'^mpnfyCs, where rrip is the mass of a proton and we have neglected an order unity 
correction arising from the non- sphericity of the flow. We refer to the surface layer of the 
disk within which photoionization heats disk gas to ~10'^K as the "base of the wind." 
At the base of the wind, a balance between photoionization and radiative recombination 
yields a number density of ~ {^i/4:7rD'^Y^'^{aj-ccd)~^^'^ , where Orec = 2.6 x lO^^'^ cm^ 
s~^ is the Case B radiative recombination coefficient for hydrogen at a temperature of 
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10^ K, and we have used the fact that the optical depth to ionizing photons is unity. 
Hence, M ~ {An /a^cc^^'^'m'pCsd^^'^i^i/ D'^Y^'^ . Though this expression is only good to 
order of magnitude, its coefficient is in fortuitously good agreement with more detailed 
wind mo dels .'1^'^ Observations yield 10^°-^ s^^ Lyman continuum photons in the central 
parsec'^^, corresponding to $j_49 = 63. Using D = Ipc, ($i = 8. The central 

concentration of S-stars within 0.01 pc, which we estimate to contribute $,,49 = 0.2 from 
each of ten approximately IOMq stars comparable to the second-most luminous Trapezium 
star,'^ yields a small total of $4,49 = 2 but in a more concentrated region. At the current 
position of the cloud, these stars contribute {^iAg/ D'^J^^'^ ~ 230 for = 6 x 10^'^pc. 
At D = 0.04pc, this number is 35. For d = lOAU and smaller, mass loss from these 
ionizing fluxes dominates over FUV-driven mass loss. Using an intermediate value of 
($i,49/^pc)^^^ = 100, ~ 3 X IO-^Mq yr-i(rfout/10AU)3/2. On the cloud's original orbit, 
M„ ~ IO-^Mq yr-i(rfout/10AU)3/2, allowing our nominal disk, which contains ~10 ^Mq 
between 1 and 8 AU, to survive for ~10^ yr. Disk masses several times larger are plausible, 
and hence a proto-planetary disk could have survived until the current time on the star's 
birth orbit in the ring. 

Dynamics of stripped gas. Currently, gas farther than ~12 AU from the star (for 
m* = O.3M0; see Figure 2) was originally ejected in the photoevaporative wind. This 
ejected material (which starts in a ring-like configuration) itself undergoes tidal stripping. 
Along the star's original orbit, the extent of the wind moving at ~ 10 km s~^ is set by the 
original 24 AU tidal radius. Since gas requires only 10 years to travel 24 AU at 10 km s~^, 
this wind scale applies even if the wind region was disrupted by close stellar encounters or 
more distant encounters with the black hole at some point in the past. 

As the star plunges toward the Galactic centre, its disk and wind are pulled off in 
shells as the tidal radius shrinks (Figure 3). The time for a parcel of wind to travel at 
10 km s~^ from 10 to 100 AU is comparable to the 70 year time to plunge to pericenter on 
the cloud's current orbit. This wind-generated cloud in turn experiences tidal disruption. 
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By its current location, the original wind cloud will have reached an extent of a few hundred 
AU. Figure 2 illustrates this extent for = O.SMq. We note that although previous close 
encounters with the black hole may have stripped the disk to smaller than its original size, 
the disk wind is regenerated over each orbit. As long as the disk size exceeds d^sc and most 
of the disk mass remains intact, our wind calculation remains valid. If the disk is stripped 
to smaller sizes, a wind will still be blown, but with reduced M^. 

Ram pressure of the ambient gas exceeds the ram pressure of the photoevaporative 
wind when Uambvl > nv^j, where Uamb is the ambient number density of gas and is the 
star's velocity along its orbit. The characteristic density of ambient gas within the central 
1.5pc is Uamb ~ lO^cm"^.!^ Models place the density at ~ 3 x 10^-6 x 10^ cm~^ at the 
cloud's current location and ~ (1^5) x 10^ cm^^ on its original orbit .'^ Along the star's 
original orbit, ~ 700 km s^^, and the ram pressure force from the disk wind roughly 
balances ram pressure with the ambient medium at the star's tidal radius. Currently, 
Vi, ~ 2300 km s~^, so ram pressure with the surrounding medium has increased by 1-2 
orders of magnitude at comparable separations from the star, and the tidally-disrupted 
photoevaporative wind is undergoing ram pressure stripping. Nevertheless, at the current 
tidal radius, the two pressures remain in rough balance and our estimates of the mass loss 
rate are therefore valid. As the cloud continues its plunge toward the super-massive black 
hole, its outer (tidally detached) extent will be shaped by ram pressure stripping. 

Figure 4a displays the inferred ionized density of the cloud as a function of radial scale. 
From the total Br7 line luminosity, the discovery paperl^ calculates an electron density of 
2.6 X 10^((i/125ylf/)~^/^, in excellent agreement with our prediction. In Figure 4a, we also 
plot the contribution to the total line luminosity as a function of radial scale. Since the 
luminosity is proportional to ri^d^, this contribution peaks at the outer edge of the disk, 
but the contribution from the extended cloud falls off slowly, as 1/d, so that about 1/5 of 
the line luminosity comes from 100 AU scales. We note that though the majority of Br7 
emission comes from the 10-20 AU scale of the tidally expanded disk, the majority of the 
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mass in the cloud is at large scales since the cloud mass is proportional to nd? oc d. At 
these large scales, full hydrodynamic simulations including ram pressure stripping and tidal 
gravity in 3D are required to match in detail the surface brightness, shape, and velocity 
width of the observed emission. One might naively expect the surface of the cloud to be 
Kelvin-Helmholtz unstable.^ However, observations of cold fronts in X-ray clusters indicate 
that gas clouds moving at a Mach number of order unity through a hot (~ keV) ambient 
medium maintain a smooth surface, probably due to "magnetic draping" .^3 

Discussion 

We predict (Figure 4b) that the total Br7 luminosity of the cloud will increase as it 
approaches pericenter. This future evolution of a debris cloud around a low-mass star on a 
Keplerian orbit is easily distinguishable from that of a pressure-confined cloud with no self 
gravity or central mass supply. We further predict that with better resolution, the specific 
intensity of the line should increase since most of the emission is coming from a smaller 
spread in velocities than is currently resolved. 

The dust in the wind, in analogy to dust in HII regions, does not reach temperature 
equilibrium with the lO'^K gas. Gillessen et al. argue that the dust continuum emission 
at about 550 K comes from small, transiently heated grains, having a total mass equal 
to ~10"^ the total gas mass of the cloud.'^ Additional colder dust may be present. This 
relatively small dust mass may result from grain growth and settling in the proto-planetary 
disk, leaving relatively few small grains available to be lofted into the photoevaporative 
wind. The neutral, cooler disk does not contribute substantially to the observed dust 
emission. 

As demonstrated above, a ~ 3 x 1O~^M0 disk can survive on the star's original orbit 
for ~3 Myr. In the Supplementary Discussion, we discuss the likely rate at which such disks 
could be scattered onto orbits comparable to the observed gas cloud. This rate primarily 
depends on the highest mass scatterers in the ring of young stars around SgrA*. It could 
increase considerably above a minimum value 0.1% if the stellar mass function extends above 
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the observed GOMq limit (as observed elsewhere in the Galaxy), if there are intermediate 
mass black hole^^of lO^-lO^M©, or if compact star clusters (such as the observed clusters 
IRS 13N or IRS 13E) act as massive scatterers. Finally, if a 0.03 Mq disk is disrupted by 
the black hole every ~5 x lO^yr, the related events deliver ~ 6 x 10~^Mq yr^^ to the inner 
0.01 pc around SgrA*, comparable to the presently inferred accretion rate onto the black 
hole.lISl 

We note that alternative explanations invoking mass loss during rare evolutionary 
phases of stars (such as a proto-planetary nebula around an emerging white dwarf or a 
debris envelope around a merged blue straggler), have a much lower probability of occurance 
than our scenario. They share the requirement for a rare kick, but involve progenitors 
which are much less abundant than low-mass stars in the ring of young stars around SgrA*. 

The existence of proto-planetary disks in galactic nuclei has important implications: it 
could lead to a fragmentation cascade to comets, asteroids, and dust around quasars,®' and 
to bright flares due to the tidal disruption of planets,'^!' as well as transits of hypervelocity 
stars. In the Milky Way centre, the debris of proto-planetary disks offers a new probe 
of the low mass end of the mass distribution of stars which are too faint to be detected 
otherwise. 
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Figure 1: Proto-planetary disk plunging toward the galactic centre. Schematic 
illustration of our scenario. A young, low-mass star hosting a proto-planetary disk is dynam- 
ically dislodged from its original orbit (dashed circle) in the Galactic centre's young stellar 
ring (gray). As it plunges toward the supermassive black hole (SMBH) at the Galactic cen- 
tre (orbit indicated by dashed ellipse), the disk is photoevaporated and tidally disrupted, 
generating an extended dusty cloud around the star (blue circle). A trail of dust and gas is 
deposited along the star's orbit (extended blue trail). 
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Figure 2: Expansion of the infalling gas. As the star plunges toward the Galactic centre, 
its disk and the surrounding photoevaporative wind are tidally disrupted. For a fiducial stellar 
mass of m=i, = 0.3Mq, we plot the distance, d, of gas from the host star at the current epoch (solid 
line) and at the time of pericenter passage (dashed line) as a function of separation from the host 
star at the time of tidal decoupling. The dotted line represents no change. We assume that the 
material is instantaneously decoupled from the star when its separation equals the tidal radius and 
that it subsequently moves as a test particle in the gravitational field of the supermassive black 
hole. At the time of its decoupling, the test particle is started with the orbital parameters that 
the cloud had at time d/u* prior to decoupling, and the full separation between the test particle 
and the star is thereafter calculated as a function of time. Red portions of the curve represent gas 
initially in the disk, while blue portions represent gas initially in the photoevaporative cloud. For 
gas currently observed in the cloud, the tidal evolution depicted here dominates the wind structure 
at large scales, while the wind outflow at ~ 10 km s~^ dominates the structure of recently ejected 
gas near the disk rim. At pericenter, tidal evolution determines the structure of the entire wind. 
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Figure 3: Disk and wind structure. Schematic diagram of the disk and wind structure 
on the star's original orbit in the young stellar ring (a) and at the current epoch (b). On its 
original orbit, the star (black circle) hosts a protoplanetary disk (gray), which is limited in 
extent (vertical black line) to approximately one third of the star's tidal radius, dt (dashed 
circle). Photoionization efficiently launches a wind beyond desc, where gas is heated enough 
that its thermal velocity is comparable to the stellar escape velocity (vertical red line). From 
the wind launching region (blue), gas flows outward (red arrows) until it passes the tidal 
radius and is stripped from the star. At the current epoch, dt has shrunk, causing both the 
disk and the wind to tidally expand. 
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Figure 4: Observational predictions. Ionized gas in the disk wind generates Br7 line luminosities in 
good agreement with observations, (a), left axis: The cloud's ionized gas number density is inferred from 
the observed Br7 line luminosity to be'^ 2.6 x 10^(d/125AU)^'^/^ cni^'^ for emission dominated by gas at a 
distance d from the center of the cloud (blue, solid). We approximate the neutral number density in the 
disk as n — mp^T,Q{di/do)^^h^^{di/d)^, where di is the initial separation of the gas from the host star and 
we use a disk scale-height ^ O.ld. We show the ionized density, at the current epoch (black, solid) 
and at the time of pericenter passage (red, solid) for = O.SM©. Where wind gas dominates, we estimate 
71+ = Mwi'^TTiTipd'^Vw)^^ with = 3 X lO^^((iout/lOAU)'^/^M0/yr, assuming that the wind geometry is not 
yet spherically symmetric. We choose dout to be the current outer limit of the tidally disrupting disk. This 
estimate is not good at the largest separations, where tidal evolution and ram pressure stripping modify 
the gas density (see text). Where the disk itself is present, the neutral gas number density in the disk 
(gray dotted, current; orange dotted, pericenter) is 5-8 orders of magnitude larger than the ionized number 
density at the disk surface. We therefore set n+ — jii,, the number density at the base of the ionized wind 
(see text). Near the outer edge of the disk, the observed number density is in excellent agreement with 
the modeled density, given the approximate nature of this model, (a), right axis: Br7 luminosity as a 
function of distance, d, from the infalling star, calculated using Lbv-i = 2.35 x 10^^^n^(4/3)7rd'^ (in c.g.s. 
units) now (black, dashed) and at pericenter (red, dashed) for the same model. The currently measured 
line luminosity is included for reference (blue, dashed), and is in good agreement with the model. Most of 
the total integrated line emission is coming from near the current disk edge at 10-20 AU, but ^1/5 of the 
emission comes from ^100 AU scales, (b) Total Br7 line luminosity integrated over our model disk plus 
wind as a function of time, calculated as Lbt-i = 2.35 x 10^^^ x 27r J n^rf^ dd (in c.g.s. units; using 2tt to 
adjust for a non-spherical wind). We predict that the total line luminosity will increase by approximately a 
factor of 5 from its reported value (blue asterisk) as the star approaches pericenter in 2013. 
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Supplementary Figure SI: Dynamical properties of stars. Semi-major axes, peri- 
center distances, eccentricities, and inclinations for perturbing stars (m* = GOMq; a) and 
low-mass stars (test particles; b) as functions of time for our first simulation, described in 
the text. The perturbers maintain properties consistent with the young stellar ring, while a 
subset of low mass stars are delivered onto low angular momentum orbits. 
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Supplementary Figure S2: Results from N-body simulations, (a) Histogram of 
angular momenta of low mass stars at time 6 x 10^ years for the same simulation shown in 
Figure 5. The number of stars is well matched by 1-D angular momentum diffusion (equation 
S5; solid). The expectation for a relaxed system, = A^*(J/Jc)^A (dashed), is plotted for 
reference. For the observed cloud, logi^i^J / Jc) = —0.6. (b) Orbits of simulated stars (solid) 
that reach pericenter distances smaller than that of the observed cloud (dashed) in less than 
6 X 10^ years. The inner edge of the young stellar ring is plotted for reference (dot dashed). 
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Supplementary Figure S3: Angular momenta of low mass stars. We show his- 
tograms at times (a), 6x lOV (b), 2.5xlOV (c), and 10* yr (d) in a simulation containing 
24 perturbers, each with m^, = lOOM© (see text). The low- mass stellar population, which 
starts with evenly distributed angular momenta in the inner ring, is first matched by the ID 
diffusion equation (Eq. 5; solid) and then by = N^{J/Jc)'^A (dashed). 
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Supplementary Discussion 



Scattering Probability. Proto-planetary disks hosted by young stars could be stripped 
by several possible mechanisms. In the main text, we have shown that on their initial orbits 
in the young stellar ring, such disks are not depleted by mass loss in disk winds over their 
several million year lifetimes. However, the same dynamical processes that could send a 
young star plunging toward the Galactic centre could themselves lead to disk stripping. In 
particular, scattering events can strip disks and, if the star's orbital evolution is too gradual, 
earlier encounters with the central black hole could also produce stripping. To arrive on 
its current plunging orbit without losing its proto-planetary disk, a young star must: (i) 
experience sufficient orbital evolution over its lifetime due to scattering interactions at all 
distances to be placed on a highly eccentric orbit while avoiding single encounters that 
strip the disk, and then (ii) experience a final, strong encounter that substantially alters 
the star's pericenter distance, again without stripping its disk. Requirement (ii) is needed 
so that the disk's extent was not already truncated by tidal stripping at a comparable 
pericenter location during a previous plunge toward the black hole. 

We note that disk winds can be driven by photoionization even if the disk is truncated 
at closer distances to the host star than those estimated above, so requirement (ii) is not 
strictly needed for our model. However, if requirement (ii) is not satisfied, the mass loss 
rate in the disk's wind and hence the observed cloud density will be lower than estimated 
in the main text. 

In order to investigate these requirements, we first consider the stripping potential of a 
single scattering encounter. We demonstrate that a final, strong encounter with a high mass 
star or intermediate mass black hol(Pl can change the pericenter distance of a plunging 
star without stripping its proto-planetary disk. We then show that, given choices about 
the distribution of massive perturbers that are consistent with observations, evolution of a 
low-mass star to a highly eccentric orbit followed by a strong final encounter could happen 
with reasonably high probability, all without disk stripping. 
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A single scattering encounter with a perturber of mass rripert at impact parameter b and 
with relative velocity at infinity fenc changes a star's velocity by Af* ~ (G'mpert/&^)(26/fenc)- 
A given change in velocity Av^, hence requires an impact parameter b ~ 2Gmpert/(fcnc^^*)- 
Such an encounter tidally strips disk gas at radii where the induced relative velocity 
between the star and its orbiting material, At>tid ~ {Gmpcrt/b'^){2rd/b){2b/venc) exceeds the 
star's escape velocity fesc,* = (SGm^/r^)-^/^. Note that because we will be interested in 
Av^: < Vcnc, we can ignore gravitational focusing in our estimate of Aftid (encounters are not 
parabolic, as is often assumed in calculations of disk stripping during close stellar passages). 
During an encounter that changes the stellar velocity by Av^, the disk is stripped at radii 
rd > {b/2){vcsc,*/Av^), or in other words, where v^^^^^ < 2(m*/mpert)fenc(A?;*)^. We choose 
to express the location of stripping in terms of fesc,* rather than because disk winds are 
launched where the thermal velocity, fth, of gas with temperature T ~10^K is comparable 
to the escape velocity from the star. The location, r^, of this equality varies as a function of 
stellar mass. The change in stellar velocity generated by an encounter that does not strip 
gas at the disk radius with a given fesc,* is limited to 

We may now estimate the perturber mass required to change a plunging star's pericenter 
distance by of order itself without truncating the star's disk at the radius corresponding 
to fesc,* = "^^th- Since the star's orbit has an eccentricity oP e = 0.938 ~ 1, the specific 
angular momentum of its orbit around the black hole is approximately J ~ (2G'Mbh'"p)^^^, 
where Vp is the pericenter distance. A change in the star's orbital angular momentum 
of A J corresponds to a change in pericenter distance Avp ^ 2rp{AJ/J). To generate 
Avp = rp/2, a single scattering event at apocenter must change the star's apocenter velocity, 
^^apo = [G'Mbh(1 - e)/rapo]^/^, by Afapo ~ fapo/4. The cloud's current, post-scattering, 
apocenter velocity of fapo.f = (3/4)fapo = 170 km s~^ requires Af* = Afapo = "^apo.f/S. 

A single scattering encounter near the apocenter of the star's orbit around the black 
hole can therefore generate Avp = rp/2 without stripping the star's disk at the distance 
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corresponding to Vesc,* = ^th as long as the perturber has mass 
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In the numerical evaluation above, we have normalized Vth to the isothermal sound speed 
for ionized hydrogen at T = 2 x lO'^K. The exact ionized gas temperature depends on gas 
metallicity (higher metallicity leads to more effective cooling) and on the typical energies 
of ionizing photons (softer spectra generate lower temperatures). We choose two reference 
encounter velocities. First, fenc = 300 km s~^ is the difference between (4/3)fapo,f and 
the velocity of a star also reaching apocenter at r^po but having an eccentricity of 0.4, 
consistent with the average eccentricity of massive ring stars, which is observed to be'^ 
~ 0.3-0.4. Near (but not exactly at) both star's apocenters, a favorably oriented encounter 
could have this velocity and reduce the angular momentum of the star by approximately 
the amount estimated above. Second, a perturber from an isotropic population would 
typically have fenc ~ {GMBu/r ^poY^^ = 670 km s~^. From equation (S2), a member of the 
stellar population of the young ring could serve as the final perturber for a star with mass 
m^: < O.3M0, given a favorable, but not implausible, encounter geometry. A solar mass star 
likely requires an encounter with a more massive perturber, perhaps an intermediate mass 
black hole with mass Mimbh ~ IO'^Mq residing outside of the young disk. Such black holes 
may be present in the Galactic centre, and are not vilating existing dynamical constraints.'^ 
In fact, Merritt et al. (2009)'^ showed that the orbits of the S-stars could be a natural 
consequence of the interaction between the nuclear star cluster and an intermediate mass 
black hole (IMBH) with a mass > 1.5 x 10'^ Mq on a mildly eccentric orbit over a period of 
a few Myr. Again, equation (S2) is not a strict requirement for our model, but if it is not 
satisfied, M^, is reduced. 

Having established that a final strong encounter can cause a disk-laden star to plunge 
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toward the black hole without stripping the disk, we now consider the probability of such 
a history for a star in the Galactic centre. We calculate the rate of delivery of low mass 
stars onto highly eccentric orbits, the probability of a final encounter that substantially 
changes the star's pericenter distance, and the likelihood that a star could be delivered 
onto a plunging orbit without experiencing a disk-stripping encounter during its lifetime. 
We make choices about the mass function of stars in the region that are simultaneously 
optimistic for our theory and consistent with observations. 

First, what is the rate of delivery of low-mass stars onto highly eccentric orbits? 
Over its lifetime T, each low-mass star is scattered multiple times by ring stars of mass 
mpert at impact parameters b > fomin, where 6min is the typical minimum encounter 
distance after an elapsed time T. The value of 6min is set by the encounter rate 
/pcrt(^) = [^pcrt/(C^-R^^)]^^^'^rei, such that /pert 7" = 1 for 6 = bmin- To avoid disk disruption, 
6min must be large enough that disk disrupting encounters are not common. We quantify 
this constraint below. Here, A^pert is the number of perturbing ring stars, Vrei is the typical 
relative velocity between perturbers and the low-mass star, R = rapo = 8400AU is the 
local ring radius, h ~ cr^,/Q is the ring scale height (with h/R 10°),'^ is the velocity 
dispersion of ring stars, and = {GMbu/ R^Y^"^ is the angular orbital velocity of the ring. 
We include an order unity factor, (, to more accurately approximate the number density of 
perturbers. We will show below that C ~ 5 given our fiducial choices for A^pert and R. The 
measured velocity dispersion cxj, = hQ = 120 km s~^ corresponds to expected eccentricities 
of V2aJ{Rn) = 0.25, in rough agreement with observed eccentricities of e* = 0.3 — 0.4. 
As a compromise between these measured values, we define a stellar random velocity 
■^^rand = V^cT* ~ e^,i?fi = 200 km s~^, and we set the typical relative rencounter velocity 
Vrei = Vrand duriug the low-mass star's early orbital evolution, where we have plugged in 
6=1, = 0.3. Once the star has reached a highly eccentric orbit, v^ei ~ ^R — fapo > cr* near 
apocenter, and we adopt Vj-d = ^IR. 

Diffusion of angular momentum. Transport of stars onto low angular momentum 
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orbits proceeds by angular momentum diffusion, as in standard loss cone calculations.'^ At 
early times, the impact parameter, ^strong, required to generate a "strong" scattering that 
alters the star's random velocity by of order itself is ^strong = ^'mpcrt/'^rand' each such 
scattering generates a fractional angular momentum change of order Vj-and/i^R)- More 
generally, as long as 6 > 6min > ^strong, each star's specific angular momentum, J ^ -Rfapo, 
is altered by of order itself on a timescale T^is = J^/D, where the diffusion coefficient 



where Jc = R^VL is the specific angular momentum of a star on a circular orbit at the 
ring radius -R, and InA = ln(0.08 pc/6min) is the Coulomb logarithm due to multiple 
gentle scatterings. Though the diffusion time for stars with angular momentum J ^ Jc 
is shorter than the relaxation time Trei = Jl/D, hy a, factor'^ of {J/JcY, the fraction of 
phase space occupied by low angular momentum orbits is also proportional to {J/Jc)"^, so 
that a population of N^, low-mass stars in the the young stellar ring produces low angular 
momentum stars a ratel^ 



Jrel 

~ fiAr.^lnA^^-^*V(^^)' 



C V ^BH / (^*VTel 

as long as the system is sufficiently relaxed. We return to the question of relaxation below. 
At early times, frei ~ ^^rand, and at late times frei ~ R^- However, we note that at late 
times, the number of perturbers available with comparable relative velocities could easily be 
a factor of RVL/v^^^-^^ ~ 3 higher given contributions from stellar remnants in the isotropic 
population, so in estimating our rates, we employ v^ei = v^anA- We note that though it 
plunges interior to the ring's extent as it approaches pericenter, the cloud nevertheless spends 
more than half of its orbital time within the ring (near apocenter). As demonstrated above, 
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disk survival requires At>tid < Vth, which translates to 6min ^ [SC^Tnpcrt'm* / {vencVth)Y^'^ ■ For 
rripert = 6OM0, Wenc = 200 km s~^, and = O.SM©, disk survival requires 6min > lOAU. At 
the limit 6min = 10 AU, In A ^ 7. Since / oc A^pert'^perf high-mass stars dominate the stellar 
orbital evolution. 

Mass function of scatterers. For the top-heavy mass function of stars inferred in the 
young stellar r'mg,^ the number of stars A^ having mass m follows dN/ dm oc 77i-0-45±o.3 
in the mass range 7Mq < m < 6OM0, and the presence of higher-mass stars is not 
excluded. Most of the ring's mass is associated with the highest mass perturbers, since 
rn?dN/dm oc w}'^^ . The inferred total mass in the young stellar ring o^^^'^ ~ IO^Mq could 
accommodate A^ ~ 170 perturbers with cl IXIRSS 'yy^pert 

= 6OM0 each. The surface density of 
observed 0/WR stars scales radially as^'' ocr~^'^, extending from 0.8" = 0.03pc= 6700AU 
to 12" = 0.5pc= lO^AU (using a distance to the Galactic centre of Rq = 8.33kpc). 
Approximately (1/10)*^'^ ~ 1/4 of the ring's 0/WR stars reside at the inner scale of the 
ring near Tapo = l" ~ 0.04pc, the location of the proposed initial stellar orbit. Hence, 
we adopt fiducial values of mpert = QOMq and A^pert = 40 in the vicinity of the cloud's 
Tapo- Integrating over the full disk (with surface density Spert (?")), we find that A'pert = 40 
perturbers lie between Ri = 6700AU and R2 = 2 x 10'^ AU. The volumetric number 
density of perturbers at i? = 8400AU is approximately Spert/(2/i) = Npert/iC'^R'^h) with 
C= (20/3)(i?o-6-/?o-6)//?o-6 = 5. 

In contrast, the mass function of observed stars both interior and exterior to the 
young stellar ring is consistent with a Salpeter/Kroupa initial mass function (IMF).— ^ We 
suppose for the sake of argument that the IMF in the young ring changes its slope to the 
Salpeter value, dN/dm ~ m~^'^^ below a stellar mass m < rrihreak = 7Mq and that this 
mass function extends down to stars of O.3M0, consistent with the observed turnover in the 
IMF in other galactic environments'^ (but see alsd^Sl). Given these assumptions, A^ oc m°'^^ 
for m > rrihreak and A^ oc m~^'^^ for m < rrihreak- Then, the inner ring would host ~12 stars 
of m = 7Mq, N = 170 stars of m = IMq, and A^ = 900 stars of m = O.3M0. The presence 
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of a few thousand young low-mass stars in the ring is consistent with hmits from X-ray 
observations.'^ These choices generate dehvery rates of / = 6 x 10^^ yr^^ for = Mq and 
3 X 10^^ yr^^ for m^, = O.SMq. Scattering on binary stars further enhances our calculated 
rate. 

N-body simulations. However, we must worry about the rate calculation in equation 
(S4) because the young stellar ring is not relaxed.l^ The angular momentum relaxation 
time Trei ~ 3 x lO^yr, and we are interested in timescales more than an order of magnitude 
shorter than this. We perform a set of N-body simulations using the Mercury6 hybrid 
orbital integrator®^ with a timestep of 10 years and and accuracy parameter of 10~^^. 
Supplementary Figures SI and S2 show results from a simulation with 40 fully massive 
perturbers, each with m^, = 60Mq, and 900 test particles. Both populations are distributed 
in semi-major axis from 6700AU to 20100AU with surface densities in semi-major axis 
proportional to r^^'^, consistent with the entire observed disk (not just the simulated 
inner part) having mass IO^Mq. Each particle starts with an eccentricity and inclination 
randomly drawn from uniform distributions between and 0.3 and between and 10°, 
respectively. The orbital nodes, arguments of pericenter, and mean anomalies are drawn 
from a uniform distribution between and 27r. In this simulation, 6 test particles (low mass 
stars) are delivered to orbits with pericenters less than or equal to that of the infalling cloud 
between 4 and 6 Myr. The massive perturbers are only mildly excited from their initial 
disk configuration. In additional simulations (not shown) starting with lower eccentricities 
and inclinations, the perturbers self-excite to values comparable to those observed on our 
timescale of interest. The number of stars delivered before 6 Myr onto orbits with pericenter 
distances less than Vp = 270AU varies substantially from simulation to simulation and can 
include 0. To elucidate the orbital evolution of the stars, we perform a longer simulation, 
this time with 24 perturbers having m^, = IOOM0 and 940 test particles. The results are 
presented in Supplementary Figure S3. At early times, the distribution of angular momenta 
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is well matched by a 1-D diffusion equation, with the number plotted 

N (^logio = (47rDt)-'/'e-(^-^=)'/(^^*) Jln(10)iV,A , (55) 

where A is the binsize. To calculate D, we use equation (S3) with R = 10^ AU, h/R = 10°, 
C — and InA = 7. At late times, phase space fills and N — N^{J/ Jc)^A, as expected. 
At 6 Myr, the transition between these regimes is marginal, and small number statistics 
generate a large variation in outcomes given our fiducial parameters. 

The expected closest approach distance, b„im, and the rate of orbital evolution, 
/, are related by / ~ r-V-^*(6strong/&min)^ In A, where p = v^^^J{QRf = e^. 
Rearranging, 6min ~ 6strong(pA''*/A'de/)^^^(ln A)^/^, where Ndei = fT is approximately 
the number of stars delivered onto plunging orbits over the age of the system, T. 
For simplicity, we set InA = 7 since variations in b„im considered here will not 
substantially alter this value. These choices give bmin ~ ^strong if all stars are delivered 
onto plunging orbits over the system age T and femm is larger if only a fraction 
of the population is delivered. Disk disruption is avoided if the delivery fraction 
Ndei/N, < {l/8){lnA){mpe^/m,)vy{vre.nd^^R^) = 0.01(mpert/60Mo)(0.3Mo/m,), where 
we have chosen to use the expression for early times since low velocity encounters arc more 
disruptive. For reference, in the simulation shown in Supplementary Figure S2, 2/6 of 
the stars delivered onto plunging orbits at ages between 4 and 6 Myr do not experience 
disruptive encounters. 

Hence, A^*/300 = 0.5 solar-mass stars and about A^*/100 = 10 stars with m* = 0.3Mq 
can be delivered over the disk lifetime T without expecting a disk disrupting encounter. For 
T — 6 Myr, the rate / above generates delivery of about A^*/50 of each stellar population, 
for a total of 3 solar mass stars and 20 0.3Mq stars. These numbers are comparable within 
the errors in our calculation, and disks can be safe from disruption. 

Other scatterers. What about disruption from encounters with lower-mass stars or with 
high mass perturbers in an isotropic population? The rate of encounters generating a fixed 
impulse Av to the disk scales as A'pert?^pert/('^*'^i'ei) for a fixed stellar mass. In a given 
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dynamical population, high mass stars dominate disk disruption, so we need not worry 
about encounters with low-mass stars in the young disk. For an isotropic population, 
(cr^t^rei)^^ is approximately 0.13 times that for the disk population, so both the star's 
orbital evolution rate / and the disk disruption rate scale as 0.13A^pert'^perf Dynamical 
limits constrain the mass of stars and stellar remnants in the galaxy's central 0.1 pc to 
be <1O^M0.151 The expected number of low-mass stars and stellar to intermediate mass 
black holes in this region is not clear .f^^ The diffusion of lighter stars out of the ring may 
have led to the unusually top-heavy mass function of stars that are found in the ring at 
present ,'^2'^^'^ and more generally, the isotropic distribution likely contains the remnants of 
older star formation episodes. For example, 10^ 1^© perturbers in an isotropic population 
would disrupt disks at ~1/10 the rate of our fiducial 60Mq perturbers. A large population 
of intermediate-mass black holes, in contrast, could generate higher disk disruption rates, 
but the number of such objects, if they exist, is unconstrained and need not be high enough 
to cause problems for our scenario. 

Finally, we return to the probability of a final strong encounter that altered the 
plunging star's pericenter distance from 2rp to Vp. The ratio of scattering events at b^i^ 
to weak scatterings is l/(lnA), so about 1/7 injected stars will come in with a strong 
(i.e., b ~ 6min) final scattering event. The rate of orbital evolution due to "ideal" strong 
scatterings with a favorable encounter velocity due to the encounter geometry is increased 
relative to the rate of typical strong scatterings by the ratio of the typical to the favorable 
encounter velocity (gravitational focusing triumphs) and is decreased by the fraction of time 
spent by a perturber with a favorable velocity. How likely strong scatterings are to occur at 
favorable encounter velocities depends on the details of the distribution of perturbers, but 
since we prefer encounters near perturbers' apocenters, where they spend the most time, 
favorable encounter geometries are not rare. 

If one IMBH having mass Mjmbh = 1.5 x IO^Mq resides in the central R = 0.04 pc of 
the galaxy, then the likelihood of a strong encounter at ^strong = GMimbh/o"^ = 3AU with 
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an isotropic ~ 670 km s~ , capable of sending a star plunging directly from a circular 
orbit, is approximately l/(47ri?^/3)7r6g|.j.Q^gCr^T = 0.01 per star for T = 6 x 10^ years. This 
corresponds to 9 injections over the disk lifetime for A^^, = 900. Stellar encounters with an 
IMBH at 3AU generate substantial disk stripping. The probability per star of a maximal 
non-stripping encounter with the IMBH approaches 1, and each such encounter can 
generate Af* ~ 150(0. 3MQ/m,,)^/^ km s^^. One or more IMBHs in the central region could 
substantially increase the rate of diffusion of low mass stars onto low angular momentum 
orbits. 

Summary. In conclusion, approximately 1/100 of the population of low-mass stars in the 
young stellar ring could have been delivered by the observed population of 60 Mq perturbers 
onto orbits comparable to that of the observed plunging cloud over the age of the young 
ring without disk stripping. A few thousand young low-mass stars could be present in the 
ring, consistent with limits from X-ray observations,'^ of which about a quarter likely lie in 
the inner ring. These constraints allow approximately 2 x 10^^ injections of undisrupted 
disks per year. We conclude that a tail population of low-mass (previously unobserved) 
stars could diffuse into orbits similar to that of the gas cloud around SgrA*. Is this delivery 
rate sufficient? A disk orbiting a O.SMq star must have a pericenter distance approximately 
2-3 times larger than that of the current cloud in order to remain undisrupted at desc- On 
such an orbit, the disk will remain bound indefinitely, and the wind will be regenerated 
over each orbital period, making the probability of observing such an object in the Galactic 
centre close to unity. Our calculations are sufficiently approximate that the infalling cloud 
may represent a long-lived, tidally truncated disk with a continually regenerated wind. If, 
on the other hand, we need a final kick that recently reduced the pericenter distance of the 
cloud by a factor of two, that kick needed to generate A J/ J = 1/4, as calculated above. 
Such final strong kicks will occur at approximately 16/7 of the delivery rate, where the 
factor of 16 comes from (J/AJ)^ and the factor of 7 comes from the Coulomb log. Over 
the 140 year period of the cloud, the probability that an infalling disk reaches the current 
cloud pericenter for the first time is ~0.1%. If the cloud is indeed on its first plunge after a 
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pericenter modification of a factor of two, we are somewhat fortunate to observe it. 

Our calculated probability is high enough that observing a plunging young star hosting 
a proto-planetary disk is plausible. We do not argue that the presence of such an object 
is likely a priori, but rather than the observed young cloud has properties that can be 
naturally explained given our scenario, and furthermore, our model generates predictions 
about the future evolution of the cloud that are not common to the pressure-confined 
cloud model proposed by the discovery team.l^ If our scenario is correct, photoevaporating 
proto-planetary disks with pericenter distances 2-3 times that of the infalling cloud should 
be common. 

We note that if our suggestion is correct and observations over the next few years 
confirm that the infalling gas cloud arises from the tidal evolution of a wind blowing off a 
circumstellar disk, more exotic sources for such a disk need not be considered. For example, 
a grazing collision between a star and a stellar mass black hole could both put the star 
onto a sufficiently eccentric orbit and generate debris, a small fraction of which might settle 
into a disk around the star.f^ However, we consider such a history less natural than our 
interpretation of the infalling cloud as a young proto-planetary disk. 

Other similar systems. Several previous observations in the region of the Galactic centre 
might also be attributable to proto-planetary disks. Clenet et al. found a flaring infrared 
source in the L' band at a projected separation of ~0.003pc = 670 AU from SgrA*.^'^ The 
emission is not likely to be stellar since no i^-band counterpart was detected (implying 
{K — L') > 3) and the source appears extended. If the emission is associated with dust, 
it corresponds to a dust temperature of ~ lO^K with a total luminosity of ~10Lq. It is 
possible that this source is of similar origin to the observed infalling cloud considered here, 
for which Kg, L', and M band fluxes imply emission fromi' dust with a temperature of 
550 ± 90K and a total luminosity of ~5L0. 

Fritz et al. found that 9 out of 15 objects detected in the Kg and if-bands in the star 
complex IRS 13E (~0.14pc = 3 x 10^ AU from SgrA*, with a core radius of about 2500AU) 
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are very red, consistent with them being warm dust clumps.'^ In principle, these dust 
clumps could be formed by the evaporation of proto-planetary disks around young stars in 
this complex. 

Finally, Muzic et al. reported another star complex IRS 13N at a similar projected 
separation from SgrA*. ™' This complex has a radius of ~2000 AU and includes extremely 
red sources with colors of either dust-embedded stars older than a few Myr or extremely 
young stars with ages < 1 Myr. The latter interpretation is supported by the fact that six 
of the sources are close in projection and show very similar proper motion whose coherence 
is not expected to survive over an orbital time. Such young stellar objects could naturally 
host proto-planetary disks of the type required to explain the infalling cloud considered 
here. 

Compact star clusters such as IRS 13N (or IRS 13E), with mass estimates of a few 
thousand solar masses, could also serve as massive perturbers able to scatter a low-mass 
star onto the orbit of the plunging cloud. 
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